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Topological susceptibility for the SU(3) Yang-Mills theory 
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We present the results of a computation of the topological susceptibility in the SU(3) Yang-Mills theory 
performed by employing the expression of the topological charge density operator suggested by Neuberger's 
fermions. In the continuum limit we find Tq\ = 0.059(3), which corresponds to \ = (191 ± 5 MeV) 4 if Fk is used 
to set the scale. Our result supports the Witten-Veneziano explanation for the large mass of the rj' . 



1. Introduction 

The topological susceptibility in the Yang- 
Mills (YM) gauge theory can be formally defined 
in Euclidean space-time as 



(1) 



X = J d A x (q(x)q(0)) , 

where the topological charge density q(x) is 



q(x) 



1 



327r 2^P' 



(2) 



Besides its interest within the pure gauge the- 
ory, x plays a crucial role in the QCD-based ex- 
planation of the large mass of the rj' meson pro- 
posed by Witten and Veneziano ( WV) |1I2| . The 
WV mechanism predicts that, at leading order 
in Ni/N c , the contribution due to the anomaly to 
the mass of the E/a(1) particle is given by |1I2I3I4| 



F>1 
2N f 



X: 



(3) 



where F v is the pion decay constant. The dis- 
covery of a fermion operator [H] that satisfies the 
Ginsparg- Wilson (GW) relation [H] triggered a 
breakthrough in the understanding of the topo- 
logical properties of the YM vacuum 7 8 4 9 10 , 
and made it possible to give a precise and unam- 
biguous implementation of the WV formula 0]. 
Indeed the naive lattice definition of the topo- 
logical susceptibility has a finite continuum limit, 
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which is universal JTOj, if the topological charge 
density, defined as suggested by GW fermions, is 
employed |4I9I10| . This yields the suggestive for- 
mula 



lim 

a -> 

V-y oo 



V 



(4) 



with Q = J2 X l( x ) = n + — n - b cm g the topo- 
logical charge, V the volume, and n+ the 
number of zero modes of D with positive (nega- 
tive) chirality in a given background. Using new 
simulation algorithms |TTj . it is now possible to 
investigate the WV scenario from first principles 
for the first time. More precisely, the aim of 
the work presented here, and fully described in 
[T^] . is to achieve an accurate and reliable de- 
termination of x in the continuum limit, which 
in turn allows a verification of the WV mecha- 
nism for the rj' mass. Several exploratory com- 
putations have already studied the susceptibility 
employing the GW definition of the topological 
charge 



2. Lattice computation 

The ensembles of gauge configurations are gen- 
erated with the standard Wilson action and peri- 
odic boundary conditions, using a combination of 
heat-bath and over-relaxation updates. More de- 
tails on the generation of the gauge configurations 
can be found in Refs. |18I19| . Tabled shows the 
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Figure 1. Comparison with other computations. 

list of simulated lattices, where the bare coupling 
constant f3 = 6/<7q, the linear size L/a in each 
direction and the number of independent config- 
urations are reported for each lattice. The topo- 
logical charge density is defined as 



q(x) = -- Tr 



j 5 D(x,x) 



(5) 



with D being the massless Neuberger-Dirac op- 
erator U3, and s an adjustable parameter in the 
range |s| < 1. For a given gauge configura- 
tion, the topological charge is computed by count- 
ing the number of zero modes of D (the inter- 
ested reader should refer to for more details). 
As s is varied, D defines a one-parameter fam- 
ily of fermion discretizations, which correspond 
to the same continuum theory but with differ- 
ent discretization errors at finite lattice spacing. 
Our computation includes data sets computed for 
s = 0.4 and s = 0.0. A comparison with re- 
sults previously obtained with various implemen- 
tation of the Neuberger's operator is shown in 
Fig. H The systematic uncertainties stem from 
finite-volume effects and from the extrapolation 
needed to reach the continuum limit. The pure 
gauge theory has a mass gap, and therefore the 
topological susceptibility approaches the infinite- 
volume limit exponentially fast with L. Since the 
mass of the lightest glueball is around 1.5 GeV, 
finite- volume effects are expected to be far below 
our statistical errors as soon as L > 1 fm. In or- 
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Ai 


6.0 


12 


2452 


1.633(48) 


0.0654(22) 


A 2 


6.1791 


16 


1138 


1.589(76) 


0.0629(32) 


A 3 


5.8989 


10 


1460 


1.737(72) 


0.0696(30) 


A 4 


6.0938 


14 


1405 


1.535(63) 


0.0615(27) 


Bq 


5.8458 


12 


2918 


5.61(16) 


0.0715(22) 


Bi 


6.0 


16 


1001 


5.58(28) 


0.0707(37) 


B 2 


6.1366 


20 


963 


4.81(24) 


0.0604(32) 


B, 


5.9249 


14 


1284 


5.59(24) 


0.0708(33) 


Co 


5.8784 


16 


1109 


15.02(72) 


0.0784(39) 


Ci 


6.0 


20 


931 


12.76(95) 


0.0662(50) 


D 


6.0 


14 


1577 


3.01(12) 


0.0651(27) 


E 


5.9 


12 


1349 


2.79(12) 


0.0543(24) 


F 


5.95 


12 


1291 


1.955(79) 0.0551(24) 


G 


G.O 


12 


3586 1.489(37) 0.0596(18) 


H 


6.1 


16 


962 


2.45(13) 


0.0599(33) 


J 


6.2 


18 


1721 2.114(76) 0.0591(24) 



Table 1 

Simulation parameters and results. 
Ai-D and E-J, s = 0.4 and s = 0.0 



For lattices 
respectively. 



der to further verify that no sizeable finite- volume 
effects are present in our data, we simulated four 
lattices at j3 = 6.0 but with different linear sizes. 
The results obtained for x, reported in Tabled 
show no dependence on L within our statistical 
errors. At finite lattice spacing, \ is affected by 
discretization effects starting at 0(a 2 ), which are 
not universal, and, in our case, depend on the 
value of s chosen to define the Neuberger opera- 
tor. The values of the adimensional quantity Tq\ 
that we obtain, where r^/a is taken from |2U|. are 
reported in Table Data, displayed in Fig. [21 
as a function of a 2 /r 2 ., show sizeable 0(a 2 ) ef- 
fects for both the s = 0.4 and s = 0.0 samples. 
For j3 < 6.0, the difference between the two dis- 
cretizations is statistically significant. A detailed 
data analysis indicate a linear dependence in a 2 
of our results within our statistical errors ■ 

3. Conclusions 

A robust estimate of x m t ne continuum limit 
can be obtained by performing a combined lin- 
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Figure 2. Continuum extrapolation of the adi- 
mensional product r^x- The s = 0.0 and s = 0.4 
data sets are represented by black circles and 
white squares respectively. The dashed lines rep- 
resent the results of the combined fit described 
in the text. The filled diamond at a = is the 
extrapolated value in the continuum limit. 

ear fit of the data. This fit gives a very good 



value of 



Xdof 



when all sets are included, and is 



very stable if some points at larger values of a /tq 
are removed. In particular a combined fit of all 
points with a 2 /r§ < 0.05 gives c = 0.059(3) with 
Xdof ~ 0-73. Since ro is not directly accessible 
to experiments, we express our result in phys- 
ical units by using the lattice determination of 
r F K = 0.4146(94) in the quenched theory [21] 
and we obtain 



X= (191±5MeV) 4 , 

which has to be compared with 



F 2 { 

— \r 

6 V 



(0) 



- 2m 



A" 



(180 MeV) 4 . (7) 



Notice that Eq. © being only valid at the lead- 
ing order in a N{/N c expansion, the ambiguity 
in the conversion to physical units in the pure 
gauge theory is of the same order as the neglected 
terms. Our result supports the Witten-Veneziano 
mechanism for explaining the large rj mass. A 
measurement of the connected expectation values 
(Q 2n )c with the methods presented here would 



provide interesting informations about the depen- 
dence of the free energy density on the 9 angle, or 
equivalently on the probability distribution of the 
topological charge Pq, putting on solid ground 
the results presented in Ref. [22] ■ Unfortunately 
much higher statistics are required in order to 
highlight the deviations from a Gaussian distri- 
bution; higher momenta of the topological charge 
distribution measured on our data are all com- 
patible with zero within large statistical errors. 
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